OZI violating eight-quark interactions as a thermometer for chiral transitions 
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This work is a follow-up of our recent observation that in the SU{i) flavor limit with vanishing 
current quark masses the temperature for the chiral transition is substantially reduced by adding 
eight-quark interactions to the Nambu — Jona-Lasinio Lagrangian with (7a(1) breaking. Here we 
generalize the case to realistic light and strange quark masses and confirm our prior result. Addi- 
tionally, we demonstrate that depending on the strength of OZI violating eight-quark interactions, 
the system undergoes either a rapid crossover or a first order phase transition. The meson mass 
spectra of the low lying pseudoscalars and scalars at T = are not sensitive to the difference in 
the parameter settings that correspond to these two alternatives, except for the singlet-octet mixing 
scalar channels, mainly the o meson. 

PACS numbers: ll.lO.Wx, 11.30.Rd, 11.30.Qc 



The proposal to include eight-quark interactions to 
stabilize the ground state of the combined three flavor 
Nambu — Jona-Lasinio (NJL) [l| and 't Hooft [1| La- 
grangians has been done in where the most 

general spin zero combinations have been worked out. 
The bosonized Lagrangian and the corresponding scalar 
effective potential and gap equations have been derived 
at leading order of the expansion, embracing effects 
of the different current quark masses to„ 7^ 7^ ttLs. 

Since then a number of physical applications of the 
enlarged model has been discussed at length. Among 
them the spectra and properties of the low lying pseu- 
doscalar and scalar meson nonets [ll|, the behavior of 
the hadronic vacuum in a constant magnetic field ^2*1, 
or at finite temperature [13i] . This latter study has been 
done in the 5[/(3) flavor limit for massless quarks. In 
particular, it has been realized that the critical tempera- 
ture Tc at which transitions occur from the dynamically 
broken chiral phase to the symmetric phase decreases in 
presence of the eight-quark interactions. The part of the 
multi-quark Lagrangian that violates the Okubo-Zweig- 
lizuka (OZI) rule is mainly responsible for this effect. 

Although eight-quark interactions are not required to 
stabilize the vacuum of the SU (2) NJL model, their effect 
on the critical temperature of the chiral transition has 
been studied as well for this case The conclusion 

remained unchanged. 

The OZI violating interactions have a strong impact on 
the type of chiral transition with massless quarks. If their 
strength is relatively weak, the phase transition restoring 
chiral symmetry at finite temperatures is of second order, 
as in the standard NJL model. However, with increase of 
the strength at T = the effective potential develops a 
minimum at the origin [T^ . In this case there is the pos- 
sibility of coexistence with another minimum obtained 
by spontaneously breaking the symmetry from a certain 
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critical value of the strength of the 't Hooft interaction. 
These starting configurations with double vacua at T = 
[1^ lead to a lowering of the critical temperature as well, 
but a new interesting feature is that one deals now with 
a chiral transition of first order. 

We stress that the usual picture associated with spon- 
taneous breakdown of chiral symmetry in NJL models 
is the one related with a maximum of the potential at 
the origin, induced by the large value of the four-quark 
coupling strength. The alternative picture that we pro- 
mote is only possible due to the Ua{^) breaking terms 
and eight-quark interactions. The latter guarantee global 
stability of the vacuum and do not alter the quality of 
mesonic spectra obtained in the conventional case. 

The reported decrease in temperature is welcome in 
view of recent lattice calculations [l3| , obtained for finite 
values of the quark masses. In this case there is evidence 
that a rapid crossover occurs at T ~ 150 MeV (note that 
this value is not reached in the framework of the standard 
NJL approach). In the massless case one expects a first 
order transition for the three flavor case (T8j-[2]|. 

In this letter we study the combined impact of realis- 
tic values for the current quark masses and eight-quark 
interactions on the occurence of extrema in the thermo- 
dynamic potential of the model and their behavior as 
functions of the temperature. Beside a rapid crossover 
the model also allows for first order transitions to occur. 
The latter is a new attribute of NJL models, induced by 
the eight-quark forces. We elucidate how the strength of 
eight-quark interactions versus the four-quark coupling 
constant gauges the occurence of the different allowed so- 
lutions. The temperature dependence of meson spectra 
for the low lying pseudoscalars and scalars is calculated 
in the rapid crossover case to illustrate the result. 

We do not consider here the effect of the Polyakov 
loop [13|-[23|' which is known to increase the transition 
temperature by ^ 25 MeV [17]. In relation with the NJL 
model the role of the Polyakov loop has been investigated 
in several papers, see £.17. |2^-[23]. The present study 
can be extended likewise. 

The central object of the present analysis are the gap 
equations at finite temperature. The corresponding ex- 
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pressions at T = were derived in ^T]\ within a general- 
ized heat kernel scheme which takes into account quark 
mass differences in a symmetry preserving way at each 
order of the long wave-length expansion |2g| , and we will 
consider from now on the isospin limit, to„ = ^ 



(3/o - A„,/i) = 0, 
2Af,(3/o + 2A„,/i 



67r2' 



(1) 



6tt 



0. 



This system must be solved selfconsistently with the 
stationary phase equations 



Ghu + Au + ^ huhs + ^ hui2hl 
lb 4 



0, 



--hl + ^h,i2hl + hl) + ^hl^0. 
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(2) 

Here A^s = M"^ — M^, A; = — mi, I = u,d,s, and 
Ml denote the constituent quark masses. The factors 
/i, are given by the average li = [2Ji{M^) + Ji{Mf)]/3, 
« = 0, 1, . . ., and represent one-quark-loop integrals 
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responsive to these changes, diminishing with increasing 
strength of the gi coupling. The point is that although 
spectra are essentially not sensitive to this correlation, 
the effective potential is. We will see that finally the pat- 
tern of the finite temperature restoration of chiral sym- 
metry depends on how strong the OZl violating forces 
are. 

In Tables 1-3 the low lying pseudoscalar and scalar me- 
son characterisitics are displayed at T = 0, (stars denote 
input). One sees from Table 1 how a smaller value of 
the parameter G is accompanied with an increase of the 
eight-quark coupling gi . 

Table 1 

Parameters of the model: m„ = m.d,ms, and A are given in 
MeV. The couplings have the following units: G (GeV~^), k 
(GeV~^), (?2 (GeV~*). We also show the values of con- 
stituent quark masses Mu = Md and Ms in MeV (only the 
case of global minima). 



Sets 


rriu 


rus 


Mu 


Ms A G 


— K 


91 


92 


a 


5.8 


183 


348 


544 864 10.8 


921 


0* 


0* 


b 


5.8 


181 


345 


539 867 9.19 


902 


3000* 


-902 


c 


5.9 


186 


359 


544 851 7.03 


1001 


8000* 


-47 


d 


5.8 


181 


345 


539 867 5.00 


902 


10000* 


-902 



with the Pauli-Villars regularization kernel [2£ 
p{tA') 1 - (1 -I- tA') exp{~tA'), 



(4) 



where A is an ultra-violet cutoff (the model is not renor- 
malizable). In the approximation considered we need 
only to know 



Jo(A/^) =A^-M^ln 1 + — 



A_2 

M2 



and 



MM' 



In 1 



A^ 

A/2 



A2 



A2 + Af2 
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The model parameters are the four-quark coupling G, 
the 't Hooft interaction coupling k, the eight-quark cou- 
plings gi,g2 (gi multiplies the OZl violating combina- 
tion), the current quark masses mi, and the cutoff A. 
The stability of the effective potential is guaranteed if 
the couplings fulfill the following inequality lOj] 



91 > 0, gi + 3^2 > 0, 



91 ^16 



(7) 



We need additional input to fix parameters. The mass 
formulae of the pseudoscalar and scalar mesons, and the 
weak decay couplings J^, fx obtained in [ll| are used for 
that. The fit shows the interesting correlation between 
G and the OZl violating coupling gi: G goes down, when 
gi increases. The mass of /o(600) is the main observable 



Table 2 

The masses, weak decay constants of light pseudoscalars (in 
MeV), the singlet-octet mixing angle 6p (in degrees), and the 
quark condensates (^uuj, {ss'j expressed as usual by positive 
combinations in MeV. 







rriK 


m^i m^i' 




fx 


e. 


(u?i) — 




a 


138* 


494* 


480 958* 


92* 


118* 


-13.6 


237 


191 


b 


138* 


494* 


480 958* 


92* 


118* 


-13.6 


237 


192 


c 


138* 


494* 


477 958* 


92* 


117* 


-14.0 


235 


187 


d 


138* 


494* 


480 958* 


92* 


118* 


-13.6 


237 


192 



Table 3 

The masses of the scalar nonet (in MeV), and the correspond- 
ing singlet-octet mixing angle 6s (in degrees). 



Sets 


'7^ao(980) 


if* (800) 


"^/o(600) 


TO/o(980) 




a 


963.5 


1181 


707 


1353 


24 


b 


1024* 


1232 


605 


1378 


20 


c 


980* 


1201 


463 


1350 


24 


d 


1024* 


1232 


353 


1363 


16 



The generalization to finite temperature of these ex- 
pressions occurs in the quark loop integrals Jq, Ji. Intro- 
ducing the Matsubara frequencies 



MM') -> MM',T) = 1&t:'T 



n— — oo *■ 



(27r)3 
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and using the Poisson formula 



(8) 



(9) 



m— — OO 



where F{n) = exp[— s(2n+ l)^7r^T^], one integrates over 
the 3-momentum p , getting the result 

oo 

MM\T) = J Jp(.A2)e-*^^ 

) 

CXD / 2 \ 

l + 2i:(-l,"exp(j^) 

n— 1 ^ 



(10) 



Similarly one gets 



(11) 



It is then easy to verify that limT->cx) Jo,i{M'^ ,T) — 0, 
and at T = one recovers the starting expressions ([3]) . 

Using these expressions in Eq. ([1]), we solve the sys- 
tem (HI)-© numerically. We henceforth assume that the 
model parameters G, k, gi, 172, JTii, A do not depend on the 
temperature. As a result we obtain the temperature de- 
pendent solutions M;(T), representing the extrema of the 
thermodynamic potential. 

The result is show in Figs. [1] (set c) and [2] (set d). 
There are either one or three (M^'' = M^*\Mi*^), i = 
1,2,3 couples of solutions at fixed values of T. For set (c) 
(as well as (b)) only one branch of solutions is physical, 
i.e. positive valued. The other two have negative values 
for the light quark masses. This is in contrast with the 
SU{3) limit case with zero current quark masses, where 
one branch collapses to the origin Af„ = Md — Alg — 
for all values of the remaining model parameters and T. 

For finite values of the current quark masses the corre- 
sponding branch moves in the M„, Ms plane: only below 
a certain critical value of G does it remain positive val- 
ued at all T. This is the case shown in Fig. [2] for set (d), 
where all three branches are positive valued. 

In case (c) the physical changes in the values of the 
quark masses as functions of T are traced back to a sin- 
gle branch. One sees however that the onset of the tran- 
sition occurs at a value of T = Ta for which the other 
unphysical two branches meet and cease to exist. The 
rapid crossover occurs in the short temperature interval 
125 MeV < T < 140 MeV. 

As opposed to this scenario all three branches of set 
(d) are positive valued. Two of the branches (starting 
from the stable minimum and the saddle solution at T = 
0) merge in the physical region at a certain Tb and the 
surviving branch has a significantly lower mass value. 
Therefore in case (d) the changes involve a jump from 
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FIG. 1; Branches of M^'' (T), Afi'' (T) pairs, denoting ex- 
tremal points of the thermodynamic potential. Solid lines 
(start at T = as deepest minima), dashed lines (start as 
relative minima at T = 0) and dotted lines (saddle at T = 0) 
for the parameter set (c). Only one branch is in the physical 
(positive masses) region (solid curves). AH units are in MeV. 
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FIG. 2: The same as Fig. 1 but for set (d). All branches lie in 
the physical region and coexist up to T = T;,, from this value 
on only one branch survives, with much lower values of i\f;. 



one branch to the other, and lead to discontinuities in 
the observables. This is a first order phase transition. 

At present lattice QCD data have not unambiguously 
settled the question about the order of the chiral tran- 
sition. For physical values of the quark masses, calcula- 
tions with staggered fermions favor a smooth crossover 
transition [l^, while calculations with Wilson fermions 
predict the transition to be first order fsO*] . 

Our calculations show that OZI violating interactions 
can be important for the issue. If they are small (as in set 
(c)) the chiral symmetry restoring transition is crossover. 
However, above some critical value of the eight-quark 
coupling gi (as in set (d)) the pattern of the transition is 
switched to first order. This correlation can be used to 
set an upper (or lower) bound for the strength of the OZI 
violating coupling gi. Without eight-quark interactions 
(set (a)) the physical branch evolves as function of T 
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FIG. 3: The masses of pion, a = /o(600), r], kaon, t]' , ao, Kq 
and fg = /o(980) from bottom to top, for set (c) as functions 
of the temperature (all in MeV). 

qualitatively as in Fig. [U but the crossover takes place 
at much larger temperatures, T ~ 210 MeV, and the 
transition is much smoother. 

The masses of scalar and pseudoscalar mesons at finite 
temperature obtained for the set (c) are shown in Fig. [S] 
As can be seen there is a rapid crossover for all meson 
masses in the same temperature interval as in Fig. [1] 
Strictly speaking, neither this rapid crossover nor the first 
order transition case (d) do imply restoration of chiral or 
C/a(1) symmetry, but only the recovery of a distorted 
Wigner Weyl phase, with the minimum of the thermody- 
namic potential shifted to finite quark mass values due 
solely to flavor breaking effects. 

To get an insight in the role played by the different 
multi-quark interactions we analyze two limits, with the 
parameter set (c) as starting condition. 

Case 1: We set gi = 52 = = and remaining param- 
eters as in (c). In this limit the gap equation has only 
one solution for the considered parameter set, thus the 
system is in a distorted Wigner- Weyl phase. 

Case 2: We set k = and all other parameters fixed as 
in (c). In this case there is no J7a(1) breaking, but OZI 



violating effects are present. We verify that in this limit 
the gap equation has also only one solution, being again 
in a distorted Wigner- Weyl phase. 

Thus the spontaneous symmetry breakdown seen in 
the full set (c) at T = (and also in sets b and d) is 
driven exclusively by the 't Hooft interaction strength k. 
We wish not to include case (a) in the present discussion, 
as it violates the stability conditions of ([7]). 

In conclusion, the present study reveals that eight- 
quark interactions, whose effects are in an almost "dor- 
mant" state as far as the low lying scalar and pseu- 
doscalar mesonic spectra are concerned, are of great im- 
portance in the study of temperature effects in chiral 
multi-quark interaction lagrangians. It turns out that 
the mesonic spectra built on the spontaneously broken 
vacuum induced by the 't Hooft interaction strength, as 
opposed to the commonly considered case driven by the 
four-quark coupling, undergo a rapid crossover to the un- 
broken phase, with a slope and at a temperature which 
is regulated by the strength of the OZI violating eight- 
quark interactions. This strength can be adjusted in con- 
sonance with the four-quark coupling (keeping the re- 
maining model parameters fixed), and leaves the spectra 
unchanged, except for the sigma meson mass, which de- 
creases. This effect also explains why in the crossover re- 
gion the sigma meson mass drops slightly below the pion 
mass. A first order transition behavior is also a possible 
solution within the present approach. Additional infor- 
mation from lattice calculations and phenomenology is 
necessary to fix finally the strength of interactions. We 
expect that the role of eight-quark interactions are of 
equal importance in studies involving a dense medium 
and extensions of the model with the Polyakov loop. 
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